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^\ Abstract 



The complete next to leading logarithmic (©(q^)) corrections to the spin 



■ dependent weak deep inelastic structure functions 51, 53 and g^, are calculated 

using dimensional regularization within the HVBM method. Analysing the 
quark and gluon initiated contributions to these corrections for different values 
for the quark masses, a consistent factorization prescription for spin dependent 
quark distributions, which safely removes soft contributions, is defined. It is 
shown that within this scheme, quark initiated corrections are comparable in 
magnitude to those of gluonic origin, even though their contributions to the 
moments are small. 
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Introduction 



The very recent release of new and more precise data on polarized deep inelastic scattering 
experiments |1| have considerably increased the interest in many aspects of the internal spin 
structure of nucleons and their related experiments 0. Among them, a growing number of 
parametrizations for polarized parton distributions, designed to reproduce the experimental 
data, and cross sections for spin dependent processes, thought as further constraints for the 
distributions and cross check for the data, are available at present ^-§1]. 

However, the increasing precision of the experiments and the alternative of having large 
contributions from gluon initiated processes, require both parton distributions and cross 
sections being obtained at least at order as and with a consistent prescription for the re- 
gularization and substraction of soft contributions. It has been shown Q, in the context of 
the photon-gluon fusion process, that a careful analysis of the regularization procedures for 
collinear singularities is essential for isolating and substracting soft contributions. This ana- 
lysis not only sheds light on the controversy about this point, but emphasizes the importance 
of working in a consistent factorization scheme. The same kind of analysis, but for weak 
boson-gluon fusion processes, is clearly more involved due to the presence of different mass 
scales in the same process and must be carefully carried out. 

In recent years, several works have been presented in connection with spin dependent 
weak structure functions [0-0. Among them, those adressing to 0{as) corrections, only 
deal with gluon initiated corrections, and, in most cases, without worring about the proper 
substraction of soft contributions. In the present paper we calculate the complete next to 
leading logarithmic -quark and gluon initiated- corrections to the spin dependent weak deep 
inelastic structure functions, using dimensional regularization |T0[ and treating 75 and e^j^p^ 
according to the original proposal by t'Hooft and Veltman and systematized by Breitenlohner 



and Maison This method, hereafter referenced as HVBM, has all the advantages of 

dimensional regularization, can be implemented at higher orders of perturbation theory and 
is free from inconsistencies. In addition, it has been shown that soft contributions to different 
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processes involving massless quarks can clearly be identified and substracted in this scheme 



T^. The presence of non zero masses, eventually large compared to the factorization scale 
fifact, changes this situation requiring a different factorization prescription. This problem is 
also addressed. 

Finally, we discuss the phenomenological consequences of the alluded corrections in the 
structure functions, discriminating contributions from different origins and using parton 
distributions, defined with the above mentioned factorization prescription, which fit the 
available electromagnetic data. 

Structure functions and parton distributions 



In this section we begin establishing the definition of the weak structure functions that 
will be used throughout this paper and also our notation for parton distributions beyond 
leading order. 

The spin dependent component of the hadronic tensor, in the case of a longitudinally 
polarized target, can be written as 



F ■ q 

1 P ■ a P ■ a 

- -^^nip' - -^<f)9^ (1) 

P ■ q q^ q^ 



where q and P are the four-momentum of the exchanged virtual boson and the nucleon 
target respectively. With this definition, the structure functions have, at lowest order, the 
following form in terms of parton distributions 

9f{x)=j:cS{x)Aq,{x) (2) 

i 

where x is the usual Bjorken variable i.e. 
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Aqi are the spin dependent quark distributions, and the coefficients Cjj are given in tables 
1 and 2. The index i runs over the quark flavours and B indicates the boson exchanged in 
the process under consideration. 

The next to leading logarithmic corrections arise, at order as, from the diagrams in 
figures 1, 2 and 3. The evaluation of the diagram in figure la corresponds to the lowest 
order expressions for the structure functions (equation 2). Diagrams lb and Ic give no 
contribution for massless quarks in the Landau gauge. The interference between la and Id, 
and also the diagrams of figure 2, produce both infrared and coUinear contributions. While 
the infrared contributions cancel each other, the divergencies of collinear origin remain and 
have to be factorized in the definition of parton distributions. The diagrams of figure 3 also 
have collinear divergencies that have to be removed in the same way. In order to deal with 
the occuring divergencies we use dimensional regularization. In the diagrams considered 
above, 75 matrices are present not only due to the helicity proyectors, but also because of 
the weak interaction vertex. As we have mentioned before, we have chosen to deal with this 
object following the HVBM proposal. Matrix elements can be straightforwardly calculated 
with the program TRACER |T^, which masters most of the intricacies of the method. 

By means of the usual proyector technique, we obtain the contributions to the struc- 
ture functions from the phase space integrated matrix elements (see appendix A for calcu- 
lational details) 



X Aj 



5(1 - ^) + ^AP,,(^)(ln^ - i) + ^A/;(^) 
27r e 2tt 



Agf(-; 

z 



^AP,,(.)(ln^-i) + ^ 
27r e ZTT 



^g'C-) (4) 



+ 

The sum in this equation runs over the flavours of those quarks and antiquarks that partic- 



ipate in the vertex. APgg and AP^^ are the usual Altarelli-Parisi evolution kernels |TJ] . /i is 
the scale introduced by dimensional regularization and 

j = -^+hiA-K--iE (5) 
where e is defined through c? = 4 — 2e, being d the space-time dimension. Ag° and Ag^ are 



the bare densities of quarks and gluons respectively. The term proportional to S{1 — z) is 
the lowest order contribution arising from diagram la, while Afg and Afg are the quark and 
gluon initiated finite non logarithmic corrections. 

Defining the scale dependent next to leading logarithmic quark distributions as 



Aqi{x,Q'^ 



IX Z 



+ 



2-K 



a 
1 dz 

X Z 



AP,,{z)[\n%-\)+Af!{z) 



AP,,(z)(ln^-l) + A/f(^) 



(6) 



the structure functions are then given by 



^ ZVr ■ Jx Z '- i -y 

"1 dz 



,x 



+— E 



X 



CfAz) Aff{z)-Af!{z) Ag{-,Q 



(7) 



This definition factorizes the poles accounting for collinear singularities and the loga- 
rithmic terms associated with the Altarelli-Parisi evolution in the quark distributions. The 
terms A// and A/f , are designed to absorbe eventual soft contributions from the non log- 



arithmic corrections. Notice that in the usual MS [|T5[ scheme this substraction terms are 
not included. 

Modern parton distributions are usually defined in what is called a variable (scale de- 
pendent) flavour number scheme [0. In schemes of this kind, quarks whose masses are 
smaller than the typical energy scale fiphys are considered massless in the evolution. For 
those quarks with masses larger than the scale, there are no associated parton distributions. 
In this paper we adopt this kind of scheme, however we keep in mind the absolute mass 
hierarchy for factorization purposes. 



0{as) corrections 



In this section we calculate the soft and hard contributions corresponding to the different 
diagrams. We begin with those related to the box diagrams of figure 3. In reference [^] this 
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kind of diagram, but for an electromagnetic interaction, was evaluated using dimensional 
regularization within the HVBM method obtaining 



(2z-l)(ln l) + 2(l-z) 



(8) 



The subscript parton means that in order to obtain the structure function, the expres- 
sion must be convoluted with the appropriated bare parton distribution. Doing this and 
comparing with the electromagnetic version of equation (4), Afg can be identified with 



A/f (^) = \ 



{2z-l){\n- — ^-l) + 2(l-z) 



(9) 



It has been shown that the last term in equation (9) has a soft origin in the case of 
massless quarks, and therefore must be factorized. This means that, for massless quarks, 
A/f must be fixed accordingly 

^'n.<,,^. = 1 - ^ (10) 

However if we are dealing with quarks whose masses are smaller than the physical scale but 
are larger than the scale that defines soft and hard phenomena, they can be considered as 
massless for the evolution (active flavour) but the last term in equation (9) corresponds to 
a hard contribution and should not be factorized, i.e. 

^'^>,^... = (11) 

Of course, this implies a non vanishing gluonic contribution to the first moment of gf"^ only 
for light flavours (m < /i/act) independently of which are active (m < jjiphys)- 

The preceeding discussion defines the gluonic part A/f of our factorization prescription. 
A/f will be fixed after we evaluate the quark initiated processes. In the following we calculate 
the analogue of equation (9) for the weak structure functions and verify that our prescription 
factorizes soft contributions. 

For the finite non logarithmic contribution Afg is identical for both the electromag- 
netic and weak structure functions. The effect of the weak vertex is completely absorbed in 



the coefficient Cj^. If the exchanged boson does not produce change in the quark flavour, the 
factorization prescription works in complete analogy to the electromagnetic case. However, 
if there is change of flavour and one of the quarks has a mass larger than the factorization 
scale, whereas the other does not, then there is a soft contribution arising from the lighter 
but not from the heavier. The factorization prescription accounts for this as it should. 

For and (74, Afg vanishes. Due to equation (11), the contribution associated with the 
heavier flavours via Afs also vanishes. When the exchanged boson is a Z'', the contribution 
from the lighter quarks and antiquarks cancel between each other leading to a vanishing 
gluonic contribution to the structure functions. However, for flavour changing bosons there 
is no such cancellation. For example, if we deal with a and in a 4-flavour scheme, only 
d,u, s and c participate, and using equation (7) we have 

9ri^,Q')Uiuon.c =^/'^ Cr[-^!iz)]Ag{^,Q') 

Zn Jx z . J- L J z 

t=a,u,s 



dz Ag{-,Q') (12) 

ZTl Jx Z Z 

This result agrees, if the limit m^/Q'^ ^ is taken, with that of reference §], where the 
gluonic contribution to the structure functions was evaluated using quark masses and a 
transverse momentum cut-off as regulator. 

The quark initiated non logarithmic correction Afq comes from the processes in flgures 
1 and 2. The amplitudes associated to diagrams 2a and 2b are straightforwardly evaluated 
and integrated in the corresponding phase space. Calculated in the Landau gauge and using 
the HVBM method, the vertex correction amounts to replace 



-2 2 

^2 7 



(13) 



47r3 V / r(l - 2e) 
As in this gauge the quark self energy vanishes at order for massless quarks, there is no 
quark wave function renormalization, all the effect of the interference between diagrams la 
and Id is the replacement of the delta function in equation (4) by the delta times the factor 
between brackets in equation (13). Adding this contribution to those of diagrams 2a and 
2b, and isolating the logarithmic and pole terms, we identify the Afq term. 



For Qi the non logarithmic correction results to be 



+2 + ^-(^ + 



2 

1 - ^ / . 2 VI - zj+ 1 - ^ 
16 



3.1-^) (14) 



The last term can be traced back to have a coUinear origin for massless quarks, so this 
defines the correspondent tilde-term 

A7w.. = -y(l-^) (15) 
Analogously, we fix for heavy quarks 

^'m>,,^.. = (16) 

Equations (15) and (16) complete the definition of the factorization prescription we adopt. 



This prescription agrees with the one used in reference |T2[ for massless quarks but is clearly 
different for massive quarks. 

Finally, for and the corrections are given by 

Af!{z) = Af!iz) + ^il + z) (17) 

and 

Afliz)=Amz) + ^il-z) (18) 

The next to leading order quark initiated corrections to the structure functions g^, and 
(A/3 4 — A//) are identical to those obtained for the unpolarized structure functions Fi and 
F2, due to the same tensorial structure at partonic level |1^ . 



In a scheme where the regularization is performed keeping explicitly the masses of the 
quarks and introducing a cut off in transverse momentum, the hard contributions are ob- 
tained taking into account in each process the mass hierarchy. However, with the factoriza- 
tion prescription we propose, the non logarithmic terms Af^ and A/^ are not dependent on 
the kind of boson exchanged or the relation between the quark masses and the factorization 



scale. This kind of dependence is absorbed in the process independent A/ terms showing 
clearly the universal character of factorization. 

Numerical Results 

In this section we analyse the relevance of the correction terms we have just calculated. 
All the analysis on 0{as) corrections to the spin dependent weak structure functions avail- 
able in the literature deal only with the gluon initiated corrections. Even more, most studies 
on spin dependent electromagnetic structure functions and parton distributions disregard the 
quark initiated corrections either not including them in equation (7) (A/y = 0) or approxi- 
mating them by their effect in the moment of the structure function (A/g = —^6{l — z)). 
As the contribution to the moment is of order it is usually assumed that these correc- 
tions are small, however as the Afq are non trivial functions, nothing guarantees that the 
corrections are not comparable to the others for some x interval. Of course, one can always 
define a factorization prescription in such a way that some correction terms are absorbed in 
the parton distributions (at least for certain structure functions) provided the prescription 
is implemented consistently in other processes. This, however changes radically the inter- 
pretation of parton distributions. It is important then, to have an estimate of the relative 
weight of the corrections within the scheme employed and particularly in one that factorizes 
soft contributions in each process. 

We begin comparing numerically the quark and gluon initiated corrections to the struc- 
ture functions (second and third terms in equation (7), respectively), between each other and 
with the tree level part (first term of the same equation). In order to make the estimates, 
we need a set of spin dependent parton distributions defined at order as and within the 
factorization scheme proposed. We build the set taking for the bare densities the functional 
forms suggested by the spin dilution model and fix the free parameters of it in such a way 
that the available data on spin dependent electromagnetic structure functions is reproduced. 
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This procedure reduces considerably the number of parameters to be fixed and includes in 
the analysis other theoretical and experimental ingredients which are detailed in reference 

i- 

As the next to leading order evolution kernels for the polarized case haven't been calcu- 
lated yet, in order to evaluate data at different values of one can either use leading order 
evolution kernels or consider observables with a moderate dependence as the asymmetries at 
an average value of V^GeV'^. It has been shown [jl8| that the effect of the evolution 



in the asymmetries is small when compared to the ambiguities associated with the fitting 
procedure and the experimental errors. Figures (4) and (5) show the agreement between the 
electromagnetic data and the asymmetries calculated with the set. 

In Figures (6), (7) and (8) we show the contributions to different structure functions 
X g{[x\ X g^^ ioi) and x g^^ [x\ discriminating their different origins (naive, quark initi- 
ated correction, gluon initiated correction). The figures clearly show that even though the 
moment of the quark initiated corrections are smaller than those of the gluons, which have 
it main contribution in the very small x region, the quark initiated corrections are greater 
[x g{[x\ X (x)) or comparable {x gY~^ (x)) to the gluonic ones in most of the x interval. 
Compared to the naive contribution, the quark initiated correction can be as large as 20% 
of the former, while the ^ approximation ammount to 7%. Notice also, that there is a non 
vanishing gluon contribution to x gY^{x). In this structure function the gluon contributions 
associated with quarks u and d in the box diagram, cancel each other while there is no such 
cancellation between s and c due to the different factorization properties they have (equation 
11). 



Conclusions 



We have computed the complete 0{as) corrections to the spin dependent deep inelas- 
tic scattering weak structure functions introducing a factorization prescription that safely 
removes soft contributions taking into account the problem of mixing different mass scales 
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in a same process. Using a set of spin dependent parton distributions that reproduces the 
electromagnetic data, we have shown that the quark initiated corrections are larger or com- 
parable to those gluon initiated even though their moments are not. 
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Appendix A 

This appendix contains some calculational details related to the definitions of structure 
functions at partonic level. 

In addition to the hadronic tensor (eq.l), one can also define a partonic one which has a 
similar structure, but is written in terms of partonic structure functions. 



■ ixv pa ^IpP^ partonic . i nv , vartonic 

1 



f (P^ - ^q'^W - ^q^)gP^^"^''^ (19) 

p ■ q q^ 

The phase-space integrated matrix elements of either the quark or the gluon initiated di- 
agrams of figs.l, 2 and 3 give the quark and gluon components of the partonic tensor, 
respectively. The partonic structure functions must be convoluted with the appropriated 
bare parton distributions in order to obtain the usual hadronic ones. 

In order to isolate each partonic structure function it is customary to define the following 
proyectors 



2p-q 



1 



' 2(1 -e) 



(20) 



Pi 
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where 

Q 



This imphes 



(21) 



^partonic = J_ f A \M\l^ (22) 



47r 



where the spin dependent amplitude A |Mr^ = i |M+r^ - |M ' 



is defined in terms of 



those for partons whose polarization is parallel(+) or antiparallel(— ) to that of the target 
and is normahzed so that the tree-level diagram fig. la gives a 5(1 — z) . The n-dimensional 
phase-space dT within the HVBM method is given by 

for two outgoing particles, fc^ is the n-4 dimensional component of the momentum of one of 
the outgoing partons B. 



Appendix B 



The bare parton distributions inspired in the spin dilution model and used in this paper 
can be effectively parametrized as 

X Aq{x) = Aq x^" (1 - xf" {l + DqX + Eq x^) (24) 

where the parameters are given in Table (3). Notice that these parameters are not the 
free parameters of the modified spin dilution model to be adjusted!. The former include 
information about the unpolarized parton distributions. 
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Figure Captions 



Figure 1 a) Lowest order graph for DIS; b),c) and d) virtual gluon correction graphs to a). 
Figure 2 Real gluon emission corrections to la). 
Figure 3 Gluon contribution to DIS at order 

Figure 4 The spin-dependent proton asymmetry given by the fit compared to SMC, E-143 and 
earher EMC data. 

Figure 5 The spin-dependent neutron asymmetry given by the fit compared to E-142 data and 
the combined proton and deuteron SMC data. 

Figure 6 The the naive contribution and the quark and gluon initiated corrections to the spin- 
dependent structure function x gl{x) 

Figure 7 The same as figure 6 for xgY"^ {x) . 

Figure 8 The same as figure 6 for xg^"^ [x). 



Table Captions 

Table 1 Values for the coefficients C^^ defined in equation (2). 
Table 2 Axial and vector couplings for the Z° exchange. 
Table 3 Parameters of the spin dependent parton distributions. 
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Table 1: Values for the coefficients defined in equation (2). 





i = u,c,t 


i = d, s, h 


Cvi 

Cai 


1 — ^sin^Ow 
1 

2 


— \-\- ^sin^Ow 
1 

2 



Table 2: Axial and vector couplings for the Z° exchange. 





A 


B 


C 


D 


E 




1.052 


0.704 


4.085 


2.467 


17.34 


dv 


-1.066 


0.7685 


2.179 


-1.645 


0.4527 


u 


1.956 


1.504 


4.676 


-3.391 


2.981 


d 


1.035 


1.734 


2.905 


-3.051 


2.324 


s 


0.0932 


1.594 


8.694 


-2.149 


32.709 


c 


-2.171 


57.27 


173.6 


-20.36 


2.279 


9 


3.828 


0.5243 


2.287 


-2.747 


1.953 



Table 3: Parameters of the spin dependent parton distributions. 
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Figure l.a 



Figure l.b 




Figure l.c 



Figure l.d 




Figure 3. a 



Figure 3.b 



This figure "figl-l.png" is available in "png" format from: 
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